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Abstract 



We present results from many-body calculations for /3-stable neutron star 

matter with nucleonic and hyperonic degrees of freedom, employing the most 

recent parametrizations of the baryon-baryon interaction of the Nijmegen 

group. It is found that the only strange baryons emerging in /3-stable matter 

up to total baryonic densities of 1.2 fm -3 are S~ and A. The corresponding 

equations of state are thence used to compute properties of neutron stars such 

as the masses, moments of inertia and radii. We also study the possibility of 

forming a hyperon superfluid and discuss its implications for neutron stars. 
PACS numbers: 13.75.Ev, 21.30.-x, 21.65. +f, 26.60.+C, 97.60.Gb, 97.60.Jd 
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I. INTRODUCTION 



The physics of compact objects like neutron stars offers an intriguing interplay between 
nuclear processes and astrophysical observables. Neutron stars exhibit conditions far from 
those encountered on earth; typically, expected densities p of a neutron star interior are of 
the order of 10 3 or more times the density p d 4- 10 11 g/cm 3 at 'neutron drip', the density at 
which nuclei begin to dissolve and merge together. Thus, the determination of an equation 
of state (EoS) for dense matter is essential to calculations of neutron star properties. The 
EoS determines properties such as the mass range, the mass-radius relationship, the crust 
thickness and the cooling rate. The same EoS is also crucial in calculating the energy 
released in a supernova explosion. 

At densities near to the saturation density of nuclear matter, ( with number density 
n = 0.16 fm~ 3 ), we expect the matter to be composed of mainly neutrons, protons and 
electrons in /^-equilibrium, since neutrinos have on average a mean free path larger than the 
radius of the neutron star. The equilibrium conditions can then be summarized as 

u n = u p + « e , n p = n e , (1) 

where «j and rii refer to the chemical potential and number density in fm~ 3 of particle 
species i, respectively. At the saturation density of nuclear matter, no, the electron chemical 
potential is of the order ~ 100 MeV. Once the rest mass of the muon is exceeded, it becomes 
energetically favorable for an electron at the top of the e~ Fermi surface to decay into a 
u~. We then develop a Fermi sea of degenerate negative muons, and we have to modify the 
charge balance according to n p = n e + n M , and require that u e = u M . 

As the density increases, new hadronic degrees of freedom may appear in addition to 
neutrons and protons. One such degree of freedom is hyperons, baryons with a strangeness 
content. Contrary to terrestrial conditions where hyperons are unstable and decay into nu- 
cleons through the weak interaction, the equilibrium conditions in neutron stars can make 
the inverse process happen, so that the formation of hyperons becomes energetically favor- 
able. As soon as the chemical potential of the neutron becomes sufficiently large, energetic 
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neutrons can decay via weak strangeness non-conserving interactions into A hyperons leading 
to a A Fermi sea with Ma = Mn- However, one expects £~ to appear via 

e~ + n -> £~ + i/ e , (2) 

at lower densities than the A, even though S~ is more massive. The negatively charged 
hyperons appear in the ground state of matter when their masses equal /x e + fi n , while the 
neutral hyperon A appears when /i n equals its mass. Since the electron chemical potential 
in matter is larger than the mass difference mj- — m\ = 81.76 MeV, S~ will appear at 
lower densities than A. For matter with hyperons as well the chemical equilibrium condition 
becomes, 

MS" = A*S- = Mn + Me; 
MA = MH° = ME = Mw 

M£+ = Mp = Mn - Me- (3) 

We have omitted isobars A, see the discussion below. 

Hyperonic degrees of freedom have been considered by several authors, but mainly within 
the framework of relativistic mean field models or parametrized effective interactions 
Q, see also Balberg et al. || for a recent update. Realistic hyperon- nucleon interactions 
were employed by Schulze et al. recently, see Ref. [|J, in a many-body calculation in order to 
study where hyperons appear in neutron star matter. All these works show that hyperons 
appear at densities of the order of ~ 2uq. 

In Ref. however, one was only able to fix the density where S _ appears, since 
only a hyperon-nucleon interaction was employed. As soon as S~ appears, one needs a 
hyperon- hyperon interaction in order to estimate e.g., the self-energy of A. The aim of 
this work is thus to present results from many-body calculations of hyperonic degrees of 
freedom for /3-stable neutron star matter employing interactions which also account for 
strangeness S < —1. To achieve this goal, our many-body scheme starts with the most 
recent parametrization of the free baryon-baryon potentials for the complete baryon octet 
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as defined by Stoks and Rijken in Ref. 0. This entails a microscopic description of matter 
starting from realistic nucleon-nucleon, hyperon-nucleon and hyperon-hyperon interactions. 
In a recent work || we have developed a formalism for microscopic Brueckner-type calcula- 
tions of dense nuclear matter that includes all types of baryon-baryon interactions and allows 
to treat any asymmetry on the fractions of the different species (n,p, A, E~, S°, X + , and 
Results for various fractions of the above particles were also discussed. 
Here we extend the calculations of Ref. [§] to studies of /3-stable neutron star matter. Our 
results, together with a brief summary of the formalism discussed in Ref. |§, are presented 
in section ||. There we discuss the equation of state (EoS) and the composition of /3-stable 
matter with various baryon-baryon potentials. Based on the composition of matter we 
present also results for baryon superfluidity and discuss the possible neutron star structures. 

II. EQUATION OF STATE AND COMPOSITION OF /3-STABLE MATTER 

Our many-body scheme starts with the most recent parametrization of the free baryon- 
baryon potentials for the complete baryon octet as defined by Stoks and Rijken in Ref. 
|7]]. This potential model, which aims at describing all interaction channels with strangeness 
from S = to S = —4, is based on SU(3) extensions of the Nijmegen potential models || for 
the S = and S = —1 channels, which are fitted to the available body of experimental data 
and constrain all free parameters in the model. In our discussion we employ the interaction 
version NSC97e of Ref. 0, since this model, together with the model NSC97f of Ref. 
result in the best predicitions for hypernuclear observables ||. For a discussion of other 
interaction models, see Refs. 0,|1O[- 

A. Formalism 

With a given interaction model, the next step is to introduce effects from the nuclear 
medium. Here we will construct the so-called G-matrix, which takes into account short- 
range correlations for all strangeness sectors, and solve the equations for the single-particle 
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energies of the various baryons self-consistently. The G-matrix is formally given by 

(B 1 B 2 \ G{uo) \B 3 B 4 ) = (B 1 B 2 \ V \B 3 B 4 ) + 
£ (B 1 B 2 \V\B 5 B 6 ) 1 



B 5 B 6 



uj - e B5 - e Bfj + vr) 
x (B 5 B 6 \ G(u) \B 3 B 4 ) . (4) 



Here Bi represents all possible baryons n, p, A, £~, S°, £ + , H~ and 2° and their quantum 
numbers such as spin, isospin, strangeness, linear momenta and orbital momenta. The in- 
termediate states B 5 Bq are those which are allowed by the Pauli principle, and the energy 
variable uo is the starting energy defined by the single-particle energies of the incoming ex- 
ternal particles B 3 B 4 . The G-matrix is solved using relative and centre-of-mass coordinates, 
see e.g., Refs. [BLCT for computational details. The single-particle energies are given by 

e Bi = t Bi + u Bl + m Bi (5) 

where t Bi is the kinetic energy and m Bi the mass of baryon Bi. The single-particle potential 
u Bi is defined by 

u Bi = Re Yl (BiBjl = £ b- 3 + e Bi ) |BiBj) . (6) 

Bj<Fj 

The linear momentum of the intermediate single-particle state Bj is limited by the size of the 
Fermi surface Fj for particle species Bj. The last equation is displayed in terms of Goldstone 
diagrams in Fig. [l|. Diagram (a) represents contributions from nucleons only as hole states, 
while diagram (b) has only hyperons as holes states in case we have a finite hyperon fraction 
in /3-stable neutron star matter. The external legs represent nucleons and hyperons. 

The total non-relativistic energy density, e, and the total binding energy per baryon, S, 
can be evaluated from the baryon single-particle potentials in the following way 

rkP d 3 k fh 2 k 2 1 „A 
^ 2 ?/o WfW- B + 2 UB{k) ) (?) 



(8) 



where n is the total baryonic density. The density of a given baryon species is given by 

k 3 

n B = g^t = x B n , (9) 
where xb = ns/n is the fraction of baryon B, which is of course constrained by 

£a*=l. (10) 

B 

Detailed expressions for the single-particle energies and the G-matrices involved can be 
found in Ref. |§. In order to satisfy the equations for /3-stable matter summarized in Eq. 
(§), we need to solve Eqs. (f|) and (|) to obtain the single-particle energies of the particles 
involved at the corresponding Fermi momenta. Typically, for every total baryonic density 
n = tin + ny, the density of nucleons plus hyperons, Eqs. (|]) and (J5|) were solved for five 
nucleon fractions and five hyperons fractions and, for every nucleon and hyperon fraction, 
we computed three proton fractions and three fractions for the relevant hyperons. The set 
of equations in Eq. (Q) were then solved by interpolating between different nucleon and 
hyperon fractions. 

The many-body approach outlined above is the lowest-order Brueckner-Hartree-Fock 
(BHF) method extended to the hyperon sector. This means also that we consider only two- 
body interactions. However, it is well-known from studies of nuclear matter and neutron 
star matter with nucleonic degrees of freedom only that three-body forces are important in 
order to reproduce the saturation properties of nuclear matter, see e.g., Ref. for the 
most recent approach. In order to include such effects, we replace the contributions to the 
proton and neutron self-energies arising from intermediate nucleonic states only, see diagram 
(a) of Fig. U, with those derived from Ref. [jTT]] (hereafter APR98) where the Argonne Vis 



nucleon-nucleon interaction 121 is used with relativistic boost corrections and a fitted three- 



body interaction, model. The calculations of Ref. JTT[ represent at present perhaps the 
most sophisticated many-body approach to dense matter. In the discussions below we will 
thus present two sets of results for /3-stable matter, one where the nucleonic contributions 
to the self-energy of nucleons is derived from the baryon-baryon potential model of Stoks 
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and Rijken and one where the nucleonic contributions are replaced with the results from 
Ref. [|TTJ following the parametrization discussed in Eq. (49) of Ref. ||13|| . Replacing the 
nucleon-nucleon part of the interaction model of Ref. J7| with that from the Vis nucleon- 
nucleon interaction does not introduce large differences at the BHF level. However, 
the inclusion of three-body forces as done in Ref. |]TT| is important. Hyperonic contributions 
will however all be calculated with the baryon-baryon interaction of Stoks and Rijken @. 



B. /3-stable neutron star matter 

The above models for the pure nucleonic part combined with the hyperon contribution 
yield the composition of /3-stable matter, up to total baryonic number density n = 1.2 
fm -3 , shown in Fig. 0. The corresponding energies per baryon are shown in Fig. ^| for both 
pure nucleonic (BHF and APR98 pn-matter) and hyperonic matter (BHF and APR98 with 
hyperons) in /3-equilibrium for the same baryonic densities as in Fig. 0. 

For both types of calculations £~ appears at densities ~ 2 — 3n . Since the EoS of 
APR98 for nucleonic matter yields a stiffer EoS than the corresponding BHF calculation, 
£~ appears at n = 0.27 fm" 3 for the APR98 EoS and n = 0.35 fm" 3 for the BHF EoS. 
These results are in fair agreement with results obtained from mean field calculations, see 
e.g., Refs. |l]-Q. The introduction of hyperons leads to a considerable softening of the EoS. 
Moreover, as soon as hyperons appear, the leptons tend to disappear, totally in the APR98 
case whereas in the BHF calculation only muons disappear. For the APR98 case, positrons 
appear at higher densities, i.e., n = 1.18 fm -3 . This result is related to the fact that A 
does not appear at the densities considered here for the BHF EoS. For the APR98 EoS, A 
appears at a density n = 0.67 fm~ 3 . Recalling /x\ = \i n = \i p + /i e and that the APR98 EoS 
is stiffer due to the inclusion of three-body forces, this clearly enhances the possibility of 
creating a A with the APR98 EoS. However, the fact that A does not appear in the BHF 
calculation can also, in addition to the softer EoS, be retraced to a delicate balance between 
the nucleonic and hyperonic hole state contributions (and thereby to features of the baryon- 
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baryon interaction) to the self-energy of the baryons considered here, see diagrams (a) and 
(b) in Fig. [I]. Stated differently, the contributions from E~, proton and neutron hole states 
to the A chemical potential are not attractive enough to lower the chemical potential of the 
A so that it equals that of the neutron. Furthermore, the chemical potential of the neutron 
does not increase enough since contributions from E~ hole states to the neutron self-energy 
are attractive, see e.g., Ref. |§ for a detailed account of these aspects of the interaction 
model. 

We illustrate the role played by the two different choices for nucleonic EoS in Fig. |4] in 
terms of the chemical potentials for various baryons for matter in /3-equilibrium. We also 
note that, using the criteria in Eq. (|j), neither the E° nor E + do appear for both the BHF 
and the APR98 equations of state. This is due to the fact that none of the E°-baryon and 
E + -baryon interactions are attractive enough. A similar argument applies to 5° and H~. In 
the latter case the mass of the particle is ~ 1315 MeV and almost 200 MeV in attraction is 
needed in order to fullfil e.g., the condition fi\ = fi E o = fi n . This has also been checked by us 
Tj l in studies of the self-energy of S~ in finite nuclei, using the recipe outlined in Ref. ||T^| . 



For both light and medium heavy nuclei, S~ is unbound with the present hyperon-hyperon 
interactions, except for version NSC97f of Ref. The latter results in a weakly bound E~, 
in agreement with the recent studies of Batty et al. |l||. From the bottom panel of Fig. |] 



we see however that E° could appear at densities close to 1.2 fm -3 . Thus, for the present 
densities, which would be within the range of energies for where the interaction model has 
been fitted, the only hyperons which can appear are E~ and A. 

In summary, using the realistic EoS of Akmal et al. [|llj for the nucleonic sector and 



including hyperons through the most recent model for the baryon-baryon interaction of the 
Nijmegen group [[!]], we find through a many-body calculation for matter in /3-equilibrium 
that E~ appears at a density of n = 0.27 fm~ 3 while A appears at n = 0.67 fm~ 3 . Due to the 
formation of hyperons, the matter is deleptonized at a density of n = 0.85 fm~ 3 . Within our 
many-body approach, no other hyperons appear at densities below n = 1.2 fm~ 3 . Although 



the EoS of Akmal et al. |LT| may be viewed as the currently most realistic approach to the 



nucleonic EoS, our results have to be gauged with the uncertainty in the hyperon-hyperon 
and nucleon-hyperon interactions. Especially, if the hyperon-hyperon interactions tend to be 
more attractive, this may lead to the formation of hyperons such as the A, S°, S + , H~ and 
S° at lower densities. The hyperon-hyperon interaction and the stiffness of the nucleonic 
contribution play crucial roles in the formation of various hyperons. These results differ 
from present mean field calculations where all kinds of hyperons can appear at the 

densities considered here. 



C. Baryon superfluidity in /3-stable matter 

A generic feature of fermion systems with attractive interactions is that they may be 
superfluid in a region of the density-temperature plane. The 1 S'o wave of the nucleon- 
nucleon interaction is the best known and most investigated case in neutron stars, and the 
results indicate that one may expect a neutron superfluid in the inner crust of the star and a 
proton superfluid in the quantum liquid interior, both with energy gaps of the order of 1 MeV 
||T7^p2[. Furthermore, neutrons in the quantum liquid interior may form a superfluid due 
to the attractive 3 P2- 3 F 2 wave of the nucleon-nucleon interaction |23| . Baryon superfluidity 



has important consequences for a number of neutron star phenomena, including glitches |24[ 



and cooling |[25|| . If hyperons appear in neutron stars, they may also form superfluids if their 
interactions are sufficiently attractive. The case of A superfluidity has been investigated by 
Balberg and Barnea using parametrized effective A-A interactions. Results for A and 
S~-pairing using bare hyperon-hyperon interaction models have recently been presented 
by Takatsuka and Tamagaki p7| |. The result of both groups indicate the presence of a A 



superfluid for baryon densities in the range of 2-4no- The latter authors also suggest that 
the formation of a S~ superfluid may be more likely than A-superfluidity. Along the lines 
followed by these authors we will here present results for hyperon superfluidity within our 
model. 

The crucial quantity in determining the onset of superfluidity is the energy gap function 
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A(k). The value of this function at the Fermi surface is proportional to the critical tem- 
perature of the superfluid, and by determining A we therefore map out the region of the 
density-temperature plane where the superfluid may exist. When the l So interaction is the 
driving cause of the superfluidity, the gap function becomes isotropic and depends on the 
magnitude of k only. It can be determined by solving the BCS gap equation 

A(*) = -- r ' dk'k' 2 Vi So (hk') A [ k>) (11) 



In this equation, is the momentum-dependent single particle energy in the medium for the 
particle species in question, \i is the corresponding chemical potential, and Vis is the effective 
pairing interaction. At this point we emphasize that using parametrized effective interactions 
in the gap equation can lead to errors. The gap equation includes diagrams also found in 
the G-matrix, and one therefore needs to calculate V systematically from microscopic many- 
body theory to avoid double counting of ladder contributions. The expansion for V can be 



found in e.g. Migdal ||28|| , and to lowest order V = V, the free-space two-particle interaction. 
Higher order terms include contributions from e.g. density- and spin-density fluctuations. In 
this first exploratory calculation we will follow Ref. |27|| and use the bare hyperon-hyperon 



interaction in Eq. (|TI]). The relevant hyperon fractions and single-particle energies are taken 
from the BHF calculations described earlier in this paper. Details of the numerical solution 
of the gap equation can be found in Ref. |22| . 

Fig. H shows the energy gap A F = A(k ( F ') as a function of the total baryon density 
for £~ hyperons in /5-stable matter for the NSC97E model. Although A may appear at 
higher densities, the 1 So A-A matrix elements of the NSC97E interaction are all repulsive, 
and therefore the energy gap for A hyperons would (to lowest order) have been zero at all 
densities, i.e. these particles would not have formed a superfluid. This is it at variance 
with the results of Ref. ||26|| , however, as remarked earlier this work employs an effective, 
parametrized interaction to drive the gap equation and therefore overestimates the A energy 
gap. Our E~ results are comparable to those of Ref. [27] which were obtained with a gaussian 
soft core parametrization of the bare interaction. 
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If taken at face value these results have implications for neutron star cooling. Since 
at low densities S~ is the only hyperon species that is present in our calculation, the most 
important contribution to the neutrino cooling rate at such densities comes from the reaction 
£~ — > n + e~ + I7 e . According to Ref. [^] the threshold density for this reaction to occur 
is at around 2Atiq. If the S~s are superfluid with energy gaps similar to what we found 
here, a sizeable reduction of the order of exp(— Ap/kT) may be expected in the reaction 
rate. If neutron stars were to cool through direct Urea processes, their surface temperatures 
would be barely detectable within less than 100 yr of the star's birth. This is at askance 
with present observations. Thus, the formation of a hyperon superfluid will clearly suppress 
the hyperon direct Urea process and cooling will most likely proceed through less efficient 
processes and bring the results closer to experimental surface temperatures. 



D. Structure of neutron stars 

We end this section with a discussion on neutron star properties with the above equations 
of state. 

The best determined neutron star masses are found in binary pulsars and all lie in the 



range 1.35±O.O4M pOf except for the nonrelativistic pulsar PSR J1012+5307 of mass M 



(2.1 ± O.8)M ||31|| . Several X-ray binary masses have been measured of which the heaviest 



are Vela X-l with M = (1.9 ± 0.2)M Q |2j and Cygnus X-2 with M = (1.8 ± 0.4)M Q |3 
The recent discovery of high-frequency brightness oscillations in low-mass X-ray binaries 
provides a promising new method for determining masses and radii of neutron stars, see 



Ref. [34]]. The kilohertz quasi-periodic oscillations (QPO) occur in pairs and are most likely 



the orbital frequencies of accreting matter in Keplerian orbits around neutron stars of mass 



M and its beat frequency with the neutron star spin. According to Zhang et al. []35| and 
Kaaret et al. [|36| the accretion can for a few QPO's be tracked to its innermost stable orbit. 
For slowly rotating stars the resulting mass is M ~ 2.2M (kHz/i/Qpo). For example, the 
maximum frequency of 1060 Hz upper QPO observed in 4U 1820-30 gives M ~ 2.25M after 
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correcting for the neutron star rotation frequency. If the maximum QPO frequencies of 4U 
1608-52 (vqpo = 1125 Hz) and 4U 1636-536 {vqpo — 1228 Hz) also correspond to innermost 
stable orbits the corresponding masses are 2.1M and 1.9M©. These constraints give us an 
upper limit for the mass of the order of M ~ 2.2M and a lower limit M ~ 1.35M Q and 
restrict thereby severely the EoS for dense matter. 

In the following, we display the results for mass and radius using the equations of state 
discussed above. In order to obtain the radius and mass of a neutron star, we have solved the 
Tolman-Oppenheimer-Volkov equation with and without rotational corrections, following 
the approach of Hartle [37j, see also Ref. [ 13| . Our results are shown in in Figs. |6| and [7]. 



The equations of state we have used are those for 

1. /3-stable pn- matter with the parametrization of the results from Akmal et al. [|lTJ made 



in Ref. ||13|| . This EoS is rather stiff compared with the EoS obtained with hyperons, see 
Fig. H The EoS yields a maximum mass M ~ 1.9M without rotational corrections 
and M ~ 2.1M Q when rotational corrections are included. The results for the mass are 
shown in Fig. |] as functions of central density n c . They are labelled as pn-matter with 
and without rotational corrections. The corresponding mass-radius relation (without 
rotational corrections) is shown in Fig. [7[ 



2. The other EoS employed is that which combines the nucleonic part of Ref. |TT[ with 
the computed hyperon contribution. As can be seen from Fig. |6|, the softening of 
the EoS due to additional binding from hyperons leads to a reduction of the total 
mass. Without rotational corrections, we obtain a maximum mass M ~ 1.3M Q whilst 
the rotational correction increases the mass to M ~ 1.4M . The size of the reduc- 
tion, AM ~ 0.6 — O.7M , and the obtained neutron star masses due to hyperons are 
comparable to those reported by Balberg et al. ||. 

There are other features as well to be noted from Fig. |6|. The EoS with hyperons reaches 
a maximum mass at a central density n c ~ 1.2 — 1.3 fm -3 . In Fig. ^| we showed that the 
only hyperons which can appear at these densities are A and £~. If other hyperons were 

12 



to appear at higher densities, this would most likely lead to a further softening of the EoS, 
and thereby smaller neutron star masses. Furthermore, the softer EoS yields also a smaller 
moment of inertia, as seen in Fig. |8|. 

The reader should however note that our calculation of hyperon degrees freedom is based 
on a non-relativistic Brueckner-Hartree-Fock approach. Although the nucleonic part ex- 



tracted from Ref . |L1| , including three-body forces and relativistic boost corrections, is to be 
considered as a benchmark calculation for nucleonic degrees of freedom, relativistic effects 
in the hyperonic calculation could result in a stiffer EoS and thereby larger mass. However, 
relativistic mean field calculations with parameters which result in a similar composition of 
matter as shown in Fig. [| result in similar masses as those reported in Fig. |6|. In this sense, 
our results may provide a lower and upper bounds for the maximum mass. This leaves two 
natural options when compared to the observed neutron star masses. If the above heavy 
neutron stars prove erroneous by more detailed observations and only masses like those of 
binary pulsars are found, this may indicate that heavier neutron stars simply are not stable 
which in turn implies a soft EoS, or that a significant phase transition must occur already 
at a few times nuclear saturation densities. Our EoS with hyperons would fit into this case, 
although the mass without rotational corrections is on the lower side. Else, if the large 
masses from QPO's are confirmed, then the EoS for baryonic matter needs to be stiffer and 
in our case, this would rule out the presence of hyperons up to densities ~ 10n = 1.2 fm~ 3 . 

Although we have only considered the formation of hyperons in neutron stars, transitions 
to other degrees of freedom such as quark matter, kaon condensation and pion condensation 
may or may not take place in neutron star matter. We would however like to emphasize that 
the hyperon formation mechanisms is perhaps the most robust one and is likely to occur 
in the interior of a neutron star, unless the hyperon self-energies are strongly repulsive due 
to repulsive hyperon-nucleon and hyperon-hyperon interactions, a repulsion which would 
contradict present data on hypernuclei ||38|| . The EoS with hyperons yields however neutron 
star masses without rotational corrections which are even below ~ 1.4M . This means 
that our EoS with hyperons needs to be stiffer, a fact which may in turn imply that more 
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complicated many-body terms not included in our calculations, such as three-body forces 
between nucleons and hyperons and/or relativistic effects, are needed. 

III. CONCLUSIONS 

Employing the recent parametrization of the free baryon-baryon potentials for the com- 
plete baryon octet of Stoks and Rijken 0, we have performed a microscopic many-body 
calculation of the structure of /5-stable neutron star matter including hyperonic degrees of 
freedom. The potential model employed allows for the presence of only two types of hy- 
perons up to densities ten times nuclear matter saturation density. These hyperons are S~ 
and A. The interactions for strangeness S = —1, S = —2, S = —3 and 5* = —4 are not 
attractive enough to allow the formation of other hyperons. The presence of hyperons leads 
however to a considerable softening of the EoS, entailing a corresponding reduction of the 
maximum mass of the neutron star. With hyperons, we obtain maximum masses of the 
order M ~ 1.3 - 1.4M Q . 

In addition, since S~ hyperons appear already at total baryonic densities ~ n = 0.27 
fm -3 ), we have also considered the possibility of forming a hyperon superfTuid. The latter 
would in turn quench the increased emission of neutrinos due to the presence of hyperons. 
Within our many-body approach, we find that S~ forms a superfluid in the 1 5'o wave, 
whereas the A — A interaction for the same partial wave leads to a vanishing gap for the 
potential model employed here. 

We are much indebted to H. Heiselberg, H.-J. Schulze and V. G. J. Stoks for many usuful 
comments. This work has been supported by the DGICYT (Spain) Grant PB95-1249 and 
the Program SCR98-11 from the Generalitat de Catalunya. One of the authors (I.V.) wishes 
to acknowledge support from a doctoral fellowship of the Ministerio de Education y Cultura 
(Spain). 
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FIGURES 




a) b) 

FIG. 1. Goldstone diagrams for the single-particle potential u. a) represents the contribution 
from nucleons only as hole states while b) includes only hyperons as hole states. The wavy line 
represents the G-matrix. 
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FIG. 2. Particle densities in /3-stable neutron star matter as functions of the total baryonic 
density n. The upper panel represents the results obtained at the Brueckner-Hartree-Fock level 
with the potential of Stoks and Rijken pj. In the lower panel the nucleonic part of the self-energy 
of the nucleons has been replaced with the EoS of Ref. For the latter, (not shown in the 

figure) positrons appear at a density 1.18 fm -3 . 
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FIG. 3. Energy per baryon in /3-stable neutron star matter for different approaches as function 
of the total baryonic density n. See text for further details. 
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FIG. 4. Chemical potentials in /3-stable neutron star matter as functions of the total baryonic 
density n. The upper panel represents the results obtained at the Brueckner-Hartree-Fock level 
with the potential of Stoks and Rijken |7j. The lower panel includes results obtained with the EoS 
of Ref. |T1. 
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FIG. 5. Energy gap Aj? as a function of the total baryon density for £ hyperons in /3-stable 
matter. 
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FIG. 6. Total mass M for various equations of state. See text for further details. 
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FIG. 7. Mass-radius relation without rotational corrections for various various equations of state. 
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FIG. 8. The maximum moment of Inertia /, in units of M Q km 2 . Same equations of state as in 
the preceeding figure. All results are for /3-stable matter and rotational corrections have not been 
included in the total mass. 
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